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It is shown that different approaches towards the solution of the Einstein equations for a static 
sphericaUy symmetric perfect fluid with a 7- law equation of state lead to an Abel differential equation 
of the second kind. Its only integrable cases at present are flat spacetime, de Sitter solution and its 
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-- - I. INTRODUCTION 

^^ The Einstein equations for static and spherically symmetric perfect fluids have been investigated by many authors 

CN [Q , [^ , [^ . An abundance of solutions has been found when no equation of state is prescribed, since then the unknown 

1—^ functions are more than the equations. The problem becomes very difficult when the fluid's pressure p (r) is defined 

, ^ I as a known function of the fluid's density p (r), both depending only on the radial coordinate r. A realistic equation 
of state is provided by the Newtonian polytropes p = ^p^'^^^'' which have been studied for many years Q, 0. In the 

*■ , limit when k ^ 00 a. softer isothermal equation of state emerges 



X 



,__l ! P = np, (1) 

> ■ 
(^ ' usually called the 7-law, because it is traditionally written as p = (7 — 1) p. Physically realistic perfect fluid solutions 

\^ I have 1 < n < c». Important special cases include dust {n = 00), incoherent radiation (n — 3) and stiff fluid {n — 1) 

^^ . where the speed of sound equals the speed of light. 

C ' Even for such a simple linear relation, few analytic solutions have been found. When n = 00 and p ^ the 

^^ , pressure vanishes, giving the case of dust. In fact, the density also vanishes and what remains is trivial flat spacetime. 

^T ' When n = — 1 the pressure and the density are constant and the solution is equivalent to a vacuum solution with 

^^ ' a cosmological constant found by de Sitter. We show in this paper that the case n = — 5 is connected to it by a 

O Buchdahl transformation 0, R|, pj and is also soluble. The case n — —3 leads to the Einstein static universe [pj. 

^Y', It is clear that the unphysical cases (n < 1) should be studied too since they produce cosmological solutions or are 

?H . connected to the physical ones by a general transformation, valid also when no equation of state is prescribed. 

PfJ' A simple solution for a general n was hidden among the solutions found by Tolman |^. In fact, he studied the 

^ , field equations under simplifying assumptions for the metric components, without imposing an equation of state. The 

pressure and the density in some of his solutions happen to satisfy Eq. (1) when certain constants are sent to zero 

^ ^ or to infinity |^. The first who systematically investigated relation (1) was O. Klein. He rediscovered the Tolman 

j^ ■ solution first for n = 3 and later for arbitrary n llfl], uM- His approach and results, published in a not-readily 

available journal, remained unnoticed for a long time. Even in Ref. ||l| his second work is mentioned as referring 

to the polytropic equation of state, which is true only for its beginning. For a third time this Klein- Tolman (KT) 

solution was found by Misner and Zapolsky |12[| . The radial dependence of the density, however, was omitted due to 

a misprint, causing additional confusion. Nevertheless, this paper became a standard reference towards which further 

rediscoveries |13[ , |H] were directed llq] . Klein also found numerically a regular solution, starting in phase space from 

Minkowski spacetime and spiralling towards the KT solution. His work was based on a certain first order differential 

equation. This solution was also rediscovered later by numeric studies of two and three-dimensional autonomous 

dynamical systems [||, ||l^, [jl^l, |1^] and further analysed. 

The question whether the regular solution, parameterized by n, has an explicit expression has never been answered 

in a satisfactory way. The persistent closure of analytic methods on the irregular KT solution suggests that there are 

no more integrable cases except it and the few solutions mentioned above. On the other hand, similar problems have 

been completely solved explicitly. For example, static dust solutions are not possible. The dust must be non-static 

or charged. In both cases the general solution has been found [|l9|, |20[|, pl[ . The cylindrically symmetric static case 

has been solved in simple functions, using two different gauges | p2|, |23f| . The relation between them was clarified in 

Ref. [p4| . The planar case follows easily from the cylindrical one [E3l or from the spherical one |jl^ , |p5[ . 

The purpose of this paper is to derive the integrable cases in a unified manner and to elucidate the mathematical 

difficulty of the problem. We show that in its heart stands the Abel equation of the second kind, whose normal form 

is 
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(2) 



Its integrable cases depend on the shape of/ (z) and are tabulated m Ref. [^. The functions that emerge are in general 
transcendental, but simplify for special values of n. The integrable cases that we find are given hy n ~ —5, —3, — 1, oo 
and the KT solution with n outside the interval (—5.83, —0.17). 

In the following three sections three different approaches are discussed which invariably lead to Eq. (2) with 
functions / (z) possessing the same general structure but with different coefficients. In Sec. II we start from the 
well-known Tolman-Oppcnheimer-Volkoff (TOV) equation S, [p^ written in a general spherical metric. In Sec. Ill 
the starting point is a differential equation derived by Klein in curvature coordinates. In Sec. IV we utilize the 
approach of Haggag and Hajj-Boutros (HH) described in Ref. [Q for a stiff perfect fluid in isotropic coordinates. Sec. 
V is dedicated to the Buchdahl transformation which supplies the integrable case n — —5. Finally, Sec. VI contains 
some discussion and conclusions. 



II. ANALYSIS OF THE TOV EQUATION 



The metric of a static spherically symmetric spacetime reads H , pM] 



(3) 



where dJl^ is the metric on the two-sphere and v, A, R depend on r. The Einstein equations are written as 

1 _, f2R" ^ R'^\ f _x\' R' 



i?2 



\ R i?2 



R 



(4) 



1 



\R I R „ , 



R \ R 



(5) 



P 



R^ 
R 



R 



(6) 



where the prime means a derivative with respect to r and units are used with SttG = c = 1. The contracted Bianchi 
identity 



p' = -{P + P)v', 



(7) 



follows from Eqs. (4)-(6) and usually replaces Eq. (6). Thus Eqs. (4), (5) and (7) determine i^,X,R,p and p. One of 
the metric functions is redundant and can be used to fix the gauge. Different coordinate systems have been introduced 
in the literature. Curvature coordinates (called also Schwarzschild coordinates) are obtained when R = r. Isotropic 
coordinates have R = re'*'/^ while polar gaussian coordinates set e^ — 1. Other coordinates are known as well. Fixing 
the gauge and the equation of state equals the number of unknowns and equations. For the time being, we proceed 
in full generality to derive the TOV equation. Let us define the so-called mass function m (r) : 



m{r) = -R{l-e-^R'^). 



Then Eq. (4) may be written as 



which integrates to 



2m' 
R^R' 



1 f" 
m{R) = - / pR^dR. 
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(8) 



(9) 



(10) 



Passing from r to R dependence, inserting e from Eq. (8) and i^ji from Eq. (7) into Eq. (5) results in the general 
TOV equation 
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(p + p) {2m+pR^ 
2R{R-2m) 



(11) 



In the case of dust p = 0. Then Eq. (11) yields pm — 0, which combined with Eqs. (9)-(10) gives p = 0. There is 
no matter and the solution is trivial flat spacetime. In the general case when p = p (p), Eq. (9) shows that Eq. (11) 
is a differential equation for m {R). Let us now introduce the variables M — ni/R, D — ^pR^ and P ~ ^pR^. Then 
Eq. (11) becomes 
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{D + P){M + P) 
{l-2M)pp ■ 



(12) 



Specializing to the 7-law equation of state and introducing t = Ini? we obtain from Eqs. (12) and (9) the autonomous 
system 



(2M ~l)Dr = D 



[n + 5)M 



1 



-D 



(13) 



D = Mr + M. 



(14) 



The case n = is excluded because it gives vanishing density. This system was derived in polar gaussian coordinates 
by Collins g and further analysed by him in Refs. [|7|, p9| |. 

Now, let us insert Eq. (14) into Eq. (13) and define x ~ M — 1/2. The result is 
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{n + \){n + 2) + 2n {n+l){n + 2) 
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A{x), 
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Eq. (8) indicates that a; < 0. The solution of Eq. (15) determines M and consequently m as functions of R. Then D 
is given by Eq. (14) which determines respectively p{R). The pressure is given by Eq. (1), while v follows from Eq. 
(7) written as 



(P + P) VR 
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(17) 



Finally, e R is found from Eq. (5), written as 



e-^R-= 1+^^ 



1 + 2RvR 



(18) 



In order to determine A (r) one should specify R (r), the most simple choice being R — r. In polar gaussian coordinates 
Eq. (18) yields an equation for R (r) with separated variables. The same is true in isotropic coordinates. 
Eq. (15) simplifies enormously when x is constant, becoming 



M 



(n + 1) +^n ^^ _ ^ 
2n 



0. 



The choice Mi — leads to m = p = p = Q, i.e. to flat spacetime. The choice 

Ma 



(n + 1) +4n 



(19) 



(20) 



gives via Eq. (14) D — M — const and A = const, p — 2D jR? which is exactly the KT solution. At the centre 
the density and the pressure have poles and diverge. This solution does not exist when —5.83 = —3 — 2\/2 < n < 
—3 + 2\/2 = —0.17 because then X2 — M2 — 1/2 is positive or vanishes (when n = —1). In the intervals —0.17 < 71 < 
and n < —5.83 the solution exists but p and p have different signs. 

When Xr 7^ we can decrease the order of Eq. (15) by the standard change of variables Xr = —y (x): 



xyvx 



Eq. (21) falls in the class 



n + 1 
2n • 



{n + 1) (n + 2) + 271 (n + 1) (n + 2) 
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y + A{x). 



[91 (x) y + .90 {x)] yx = /2 (a;) y^ + /i (a;) y + /o (a;) 



(21) 



(22) 



There is a standard procedure for the solution of such equations |2^ . It consists of two changes of variables which 
bring them to the Abel equations of the second kind given by Eq. (2) or by 



WWQ — g{C,)w -^ \. 



E = exp ( — / — dx 
51 



(23) 



This procedure is much easier when 50 = f^^d /2 = const as in Eq. (21). Namely, we have w = yE, f (z) — /0-E//1, 
5(C) = /i//o£^and 



(24) 



z (x) — / — Edx, 
J 91 



(25) 



C(x)= f ^E^dx. 
J 51 



(26) 



Since /i is a simpler polynomial than /o we shall use Eqs. (2) and (25). The use of Eqs. (23) and (26) is more 
complicated, but does not bring additional integrable cases. 

Applied to Eq. (21) the chosen alternative of the general method yields 



w — yx 



(27) 



/ \in^ + 5n + 2)2x+(n+ 1) (n + 2)1 x^^^^dx. 
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The case n = 1 (stiff fluid) is special because a logarithmic term appears in z 



= -41n|.T 



2a;' 



(29) 
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(2a: + l)(4x+l) 
a; (8a; + 3) 



(30) 



The relation z (x) is transcendental and throws / (z) out of the tables with integrable cases present in Ref. |gg|. The 
case n — — 1 also leads to a logarithmic term but its coefficient vanishes. This case is integrable. In fact, we may go 
back directly to Eq. (13) which becomes Dt = 2D if M ^ 1/2 and yields p — ~p — const. This is the well-known 
de Sitter solution. When M = 1/2 Eq. (13) is satisfied identically and we obtain formally the KT solution (20) with 
n = —1, but it has e^^R'^ ~ which is unacceptable. 
In the generic case n 7^ ±1 and Eq. (28) integrates to 



2 (n - 1) 



[{n^ + 5n + 2) 2a; - (n - 1) (n + 2)] , 



(31) 
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(n~l)(2a; + l) (n^ + 6n + l) 2a; + (n + 1)' 



[(n2 + 5n + 2)2x-{n- 1) (n + 2)] [{n^ + 5n + 2) 2a; + (n + 1) (n + 2)] ' 



(32) 



The z-x connection in Eq. (31) is transcendental, except for special values of n, and Eq. (2) is non-integrable in 
general. Let us investigate the special cases. 



When n — —3, /i divides /o and / {z) = 8z, z = — ^ {2x + 1) x~^^^. This is an integrable case, corresponding to 
the Einstein static universe. It is discussed in more details in the next section. There is also the formal n — —3 case 
of the KT solution (20). It does not exist since x — 1/4 > 0. 

One may try to simplify Eq. (31) by nullifying the coefRcients on the right. The condition n'^ + 5n + 2 = gives 
n = ^ (—5 ± vTz). This is of no good since the radical enters the power of x. The other possibility n = — 2 looks 
more promising. Then Eq. (31) becomes x'^^'^ = —Sz/A and Eq. (32) reduces to 
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This function leads to an integrable equation when the coefficient in front of z is —3/16 which is not true here. 

There are several n which convert Eq. (31) into an algebraic equation for x of fourth order or lower. It can 
be solved explicitly for x (z) and the answer replaced in Eq. (32). Third and fourth order equations appear when 
n = ±1/5, ±1/7, —1/2, —3/5. The radical structure of / (z), however, is incompatible with the tables with integrable 
/ (z). Second order equations appear when n = =Fl/3 respectively: 



4x = -5 ± V25 + 48z, 



(34) 



6zx= 17± %/l7 + 422 



(35) 



Unfortunately, the only integrable functions with square roots include the term \/z^ + zq which is not present in the 
above relations. 

Comparison between Eqs. (31)-(32) and Eq. (20) shows how complex must be the numerical regular solution, 
which starts in x, y coordinates from flat spacetime and focuses on the KT solution, following a spiral around it. 
The innocent parameter n, introduced in Eq. (1), proliferates like cancer in the process of solution, ending with 
the intricate coefficients in z (x) and f (z). It even determines the transcendental or algebraic nature of z{x). In 
conclusion, the only integrable cases found within this approach are n = oo (trivial dust solution), n = —I (de Sitter 
solution), n — —3 (Einstein static universe) and M = const (the KT solution). 

Finally, let us discuss for comparison the case of planar symmetry, which is solvable. Going to polar gaussian 
coordinates, the metric element reads 



ds^ = e^^dt^ - dr^ - R^ [dxl + F (xa)^ dx\ 



2' 
3; ' 



(36) 



where F {X2) = sin 0:2 for spherical symmetry and F (2:2) = 1 for planar symmetry |17]| , |25| . It is possible to generalize 
the TOV equation to encompass both cases. Instead of Eq. (13) one should write 



(2M -K)Dr = D 



n -f 1 
-2K + (n + 5) Af + — — Z? 



(37) 



where X = 1 or 0, corresponding to spherical or planar symmetry respectively. In the second case Eq. (37) simplifies 

" + 1 - r2 



2MMrr = 



-M^ + bMMr + aM^, 



(38) 



where a = n+ 1/n + 6, b = n + 2/n + 5. This is the analog of Eq. (15). Proceeding in the same way we obtain again 
the Abel equation (2) with y = —Mr and 



w = yM 2" 



(39) 
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(40) 



nb'^ 



(41) 



As mentioned before, Eq. (2) with / (z) 
parametric form, reads 



az + P, where a and /3 are constants, is integrable. The solution, in 



r fi 



z^Ce' - -, (42) 



a 



w = Cae^, (43) 



T = - -, , (44) 



with C being an arbitrary constant. In Ref. |17| the problem was solved in a different way, by introducing the variable 



D = D/M. Then Eqs. (14) and (37) are equivalent to 

2Dr = D(n + 7+^^D]. (45) 

This is a Bernoulli equation and is easily solved. Further details may be found in Ref. 11^ where also a connection 
with earlier work |3^, |Q on the particular cases n — 1 and ?i = 3 is established. 

III. THE APPROACH OF KLEIN 

This approach was developed in curvature coordinates where Eqs. (4), (5) and (7) simplify to 

P=l^'''''-^i^-^'')' (46) 

p=iA'e-^ + l(l-e-^), (47) 

p' = -{p + p)i^'. (48) 

Imposing the 7-law, we can integrate Eq. (48): 

p = poe-("+i)^ (49) 

Like before, the case of dust n = c», p = gives a trivial solution with constant v and A. Supposing that p 7^ 0, let 
us introduce the variables s = Por^, £, = e(""'"^^''/s and 

X - ^e-^ ^ ^. (50) 

pr^ 

Obviously x and ^ always have the same sign and are positive when the pressure is positive. 

We shall derive an equation for x, similar to Eq. (15). In Sec. II, III and IV the main functions will be denoted by 
X although they are different, in order to simplify notation and to stress the role of Eq. (22) in the whole problem. 
Let us muhiply Eq. (46) by (n + 1) /4, Eq. (47) by -1/2 and sum. The result is 

A{sx)^-{n + 3){l-e-^)C^l-n. (51) 

Let us introduce 2t = Ins. It differs by a constant from the variable in the previous section, but this is not important 
since the final equation will be autonomous. Laying temporarily aside the special case n = —3, we can express ^ from 
Eq. (51) as 

C= -[2xr + in + 7)x + n~l]. (52) 

n + 6 

Eq. (46) may be written in the following way: 

4 {sO, -{n+l) (e^ - 1) e = (" + 1) e^- (53) 



The usage of Eqs. (50) and (52) transforms this relation into an autonomous second-order equation for x. The change 
of variables Xt = —y (x) applies again, leading to 

2 (n + 3) xyy^ ^ 2{n + l)y^ + 2{n + 11) xy - {n + 1) (3n + 1) y - 4 (n + 7) x'^ 

+ {n^ + 9^2 + lln + 11) x+{n+ if {n - 1) . (54) 



This is exactly the equation of Klein derived in Refs. |10|, O], where the notation n = 2nKi + l and x — {uki + 1) xki 
was used. He studied it by series expansion and numerically, and was the first to find its regular solution (see Fig. 

(1) from Ref. [fol). When x{t) is known Eq. (52) gives S,, and Eq. (50) determines both A and p. Then Eq. (49) 

becomes an expression for v and finally p is given by Eq. (1). Modulo coefficients, Eq. (54) is the same as Eq. (21) 

and also falls in the class (22). Therefore, the procedure described in Sec. II can be applied to bring it to the form 
of Eq. (2) with 

w = yx "+3, (55) 

z = - [(n+ ll)a; + 3n-|-l]a;"^, (56) 
4a;-(n + l)^ [{n + 7) x + n - I] z 



[2 {n + 11) a; - (n + 1) (3n + 1)] [{n + 11) x + 3n + 1] 

There are no logarithmic terms in z within this approach. As with Eq. (31), x (z) is transcendental except in special 
cases. 

Eq. (54) shows that the KT solution is given here by 

.0 - ^^. (58) 

The other root x — — ^^^ leads to ^ = and e'^" ^ e^ = 0, p = np ^ oo which is unacceptable. Now, since xq > for 
any n, we must ensure that ^o > in order to have positive e~^. Eq. (52) gives ^o = n^ + 6n + 1 and we obtain the 
same conditions for the existence of the KT solution as in the previous section. 

In the case n — —1 Eq. (48) yields p ^ Po, p = ^Poi while the sum of Eqs. (46) and (47) provides the relation 
2j/ = —A. Eq. (47) is a linear equation for e^^, its solution being 

e-A = i_^ + Po^2^ (59) 

r 3 

where TOq is a constant of integration, identified as the gravitating mass. This is precisely the Kottler solution M. It 
can be used as a regular interior solution when rriQ = 0. Then it becomes the de Sitter solution. The KT solution for 
n = —1 has xq = and does not exist. 

When n = —3 Eq. (51) determines directly x: x — 1 + xi/s, while Eq. (49) gives p — poe^'^ . Then we have from 
Eq. (50) e^'^"^^ — xi + s. Eq. (46) becomes a linear equation for e'^'^ when these results are taken into account. 
There are two possibilities. When the constant xi = 0, a; = 1, which is the formal KT solution (58) for n = —3. In 
fact, it does not exist. When xi 7^ the solution is 

/ Ti \ 1/2 



xie"' = l+Ci[l + —j , (60) 

where Ci is another integration constant. A regular solution is obtained when xi = 1 and C'l — 0. Then u — Q and 

e' = TT-^- (61) 

1 + Po^ 

The last two equations represent the metric of the Einstein static universe. The pressure and the density are constant, 

P=Po, P^ -3po- 

There are two cases when the coefficients of z (x) are simplified, namely n ~ —11 and n = —1/3. In the first case 
Eq. (57) reads 



- 320/ (z) ^ 75z + 16^/''^z-3/5 + 22 x IG-^/^zi/^ (62) 

The function 

/ (z) = ciz + caz"! + C3Z«^ (63) 

is integrable for a set of (qi, 92) but (—3/5, 1/5) is not among thcni. When n = —1/3 we have 

/(.) = -l^z + 1 fl^V^-V3_ ^ ('iy\-^f^. (64) 

■' ^ ^ 64 96 V 3 y 192 V 3 y ^ ' 

The set (—1/3, —5/3) is integrable, but only when ci = —3/16 which is not the case here. 

Finally, there are few n when Eq. (56) is an algebraic equation up to the fourth order. Only the cases n = 
—5, —2, 1 are candidates for integrability because they lead to quadratic equations. Of these, n = —5, —2 yield 
radicals resembling those in Eqs. (34) and (35) and should be rejected. More interesting is the case of stiff fluid when 

X = ^ fz^ - 24 ± z^/z^ - 48^ , (65) 

2z [(z2 - 60) a; - 4] 
^(-)-- 9^-24).-8] - (^^) 

The radicals in / (z) are of the necessary type, but its structure is too complex to figure in the tables with integrable 
cases. 

Like in Sec. II, the only integrable cases found are n = —3, — l,oo and x — xq- The unsuccessful candidates for 
explicit solutions have in general different values of n in the two approaches. When they coincide, as is the case n = I, 
the reasons for rejection are different - a logarithmic term in the TOV-CoUins approach and a complicated / (z) in 
the Klein approach. 

It is interesting to compare the main variable in this section x = xk to the variables in the TOV approach, 
specialized to curvature coordinates. Eq. (8) becomes 

M=- = -(l-e-^), (67) 

and transforms Eq. (50) into 

XK = ^^. (68) 

Thus, not only x = M — 1/2, but also the above combination satisfy Abel equations of the second kind. The function 
xk resembles D, used in the planar case. It is well-known that when D satisfies the Bernoulli Eq. (45), D~^ = M/D 
satisfies a linear equation. This fact stresses once more the conclusion that the spherical case is much more complicated 
than the planar one. 

IV. APPROACH IN ISOTROPIC COORDINATES 

One can pass from arbitrary to isotropic coordinates in Eqs.(4)-(6) by putting R = re^'"^ . Let us make also the 
change s — Inr. Then 

pr\^^-\ss-\\l-\s, (69) 

pr^e^^-\l + Kus + \s+2us, (70) 

pr^e^ =. -A,, + v,s + vl. (71) 



Let us introduce next the variable 

and impose the 7-law equation of state. An expression for Vg is obtained from Eqs. (69) and (70): 

t, n + 1/1 

Vs = 4t 

nt An \ t 



(72) 



(73) 



Next, let us combine Eqs. (70) and (71) and replace in them As and Vs from the above equations. A long, but 
straightforward computation produces an autonomous equation for t 



2ntt<j 



1 



2 ,4 



-2 (1 - n) tj + - (n^ + 5n + 2) ts - 2 (n + 1) (n + 2) t^t^ - 2 (n + 1)^ t 



{n+1) +2n 



t' 



n (n + 1)^ 



8 



(74) 



A solution of this master equation determines all characteristics of the metric and the fluid. 

When n = 1 Eq. (74) is exactly Eq. (9) from Ref. fl^. In this section we generalize the HH approach to arbitrary 
n and bring it to its logical end - the Abel equation (2). We first lower the order of the polynomial in Eq. (74) by 
setting t^ = x/A and then perform the change of variables 

{V^^^~\y. (75) 

The condition a; > should be maintained throughout the calculations. Eq. (74) acquires its final form 

2nxyyx = (n - 1) y^ + [{n + 1) (n + 2) a: - (n^ + 5n + 2)] y 



+ {x - 1) An+{n+ 1) - (n + 1) x 



(76) 



It falls in the class (22) and resembles Eqs. (21) and (54), but its coefficients are different functions of n. Proceeding 
like before, we get 



w = yx 



(77) 



(n - 1) z = [{n-l){n + 2)x + n'^ + 5n + 2]x ^ , 



(78) 



/W = 



{n -l){x- 1) 4n + (n + 1)^ - (n + 1)^ x 



[{n +l){n + 2)x- (n2 + 5n + 2)] [(n - 1) (n + 2) a; + n^ + 5n + 2] ' 



(79) 



in the generic case n ^ ±1. 

The stiff fluid case leads to a logarithmic term in z (x) which makes 



/W 



2{x- l)(x-2) 
3a; -4 '' 



(80) 



non-integrable. 

The case n = — 1 is pseudo-logarithmic since the coefficient in front of In a; vanishes. We have x — l/z and 
/ (z) = 2z — 2. This case is integrable and the solution is given by Eqs. (42)-(44) with a = —(3 — 2. 

The case n = —3 leads to a gross simplification of / [z) and is also soluble. One obtains / [z) — 2z i.e. a = 2, 
/3 = 0. 

The other candidate cases may be investigated in the same manner as in Sees. II and III. The values n — ±1/3 
lead to quadratic equations with radicals of the wrong type. The case n — —2 yields 



/W 



21 18 /4 
16^~ 16 V3 



4/3 



'1/3 



16 V 3 



8/3 



-5/3 



(81) 



This equation belongs to the class (63) but again Ci ^ —3/16 and integrability is not gained. 

Finally, let us discuss the KT solution in the HH approach. When x = const Eq. (76) becomes purely algebraic 
and has two roots: xi = 1 and 

4:71 

x,^l + -—-,. (82) 

(n + l) 

When n = —1, X2 does not exist. In fact, the requirement 2:2 > shows that n must satisfy the conditions derived in 
Sec. II. The first root leads to ti = ±1/2. In the first subcase Eqs. (72) and (73) give flat spacetimc. In the second 
subcase the following line element is obtained 

ds2 = dt^ - r-^ {dr^ + r^dn^) . (83) 

The transformation r — 1/r converts this element into the usual clement for flat spacetime. 

The second root X2 represents the KT solution in isotropic coordinates, namely e^ — r"^ , e^'^ — r"^ where 

a.=±^2t±^^2, (84) 

Un+ (n+ if 



a2 = ± =. (85) 

Y'4n +(71 + 1)^ 

When n = 1, ai = ±\/2 — 2 and a2 = ±-^72. These values were found in Rcf. uM, where the KT solution was 
discovered for a fifth time. 

The integrable cases found in the HH approach coincide with those found in the Klein or the TOV-Collins ap- 
proaches. In order to make connection with the last one, we must pass in Sec. II to isotropic coordinates. We 
have 

i?' = eV2(^l + ^^, (86) 

2.T = 2Af - 1 = -R'^e-^ = -A- (87) 

We cannot obtain Eq. (74) from Eq. (15) by replacing there just x from Eq. (87) because s = Inr, while r = Ini? = 
s + A/2. The necessary additional relations are 

-A 
""^ ~ 2^2 ' 



f 2f2 



Xrr = ^~-i-. (89) 



Therefore, the HH approach combines the TOV equation in isotropic coordinates together with the linear equation 
(18) for A. 

V. BUCHDAHL TRANSFORMATION AND THE CASE TV = -5 

This transformation was found by Buchdahl [p| and rediscovered by Glass and Goldman W . Its general formulation 
refers to static perfect fluid solutions, not necessarily satisfying an equation of state. In the case of spherical symmetry 
it states that if 

ds^ ^ e^^dt^ - e^ {dr^ + r^dO^) , (90) 

is the metric of a perfect fluid solution with pressure p and density p then there is a reciprocal solution which has 
vj) = —V, Xi, = Av + \ and 
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Pb = e->, (91) 



Pb = -e-^" {p + Gp) . (92) 

The transformation is simplest in isotropic coordinates and may be applied to the results obtained in the previous 
section. When p and p satisfy the 7-law, the same is true for ph and pb but with a different parameter: 

P6 = -(n + 6)p6. (93) 

Thus, the transformation of the parameter is 71 — > — (rt + 6). It is clear that when the starting solution is physically 
realistic (n > 1), the transformed one is unphysical, — (n + 6) < —7. However, when the starting solution is unphysical 
and n < —7, then the reciprocal one is physical. When n is in the interval —7 < n < 1 both solutions are unphysical. 
The existence of this transformation teaches that unphysical solutions should not be neglected a priori and the whole 
spectrum of n must be investigated. 

Let us apply the transformation to the integrable cases found in Sec. IV. The KT solution is self-reciprocal, i.e. 
transforms into itself. The reason is that 4n-|- (n + 1) is invariant under the transformation. Taking ai^ and q;2- as 
the powers of r in the starting solution, it is easy to show that a2+ = —(^2- and ai+ ~ ai_ + 2a2- when the plus 
solution is Buchdahl transformed. Solutions with n outside the interval (—5.83, —0.17) transform between themselves. 

The case n = — 3 is also self-reciprocal, as was noticed already by Buchdahl. This explains why v = when the 
element of the Einstein static universe is written in isotropic coordinates [p2| . This is necessary to ensure the equality 
between the starting and the transformed solution. 

The interesting case is the de Sitter solution (n = —1) which transforms into an explicit solution with rt = — 5. In 
isotropic coordinates the de Sitter solution reads |33] 



e^-' = 


/l + cr^y 


e^ = 


d-cr^y'.. 



(94) 
(95) 



and p = 12/ c, p = — 12/c where c is some constant. The transformed solution has 



e- = [\^ I . (96) 



1 + cr 



4 



{l-cr-^f 



cr 
I -cr'^V 12 



4 



1 + cr-^ I c 



and, of course, pb ~ —'opb- Eq- (72) supplies the corresponding t: 



1 — CT 

t= ','2y (99) 

2 (1 + cr^) 



1 - cM 
2 (1 -H 10cr2 + c'^r'^) ' 



^b - 0/1 , in„„2 , „2^4V (100) 



From Eq. (75) we get x = At^ and y = —Artj.. When tb is plugged in these relations we obtain the parametric solution 
X (r), y (r) of Eq. (76) for n = —5: 

o ft 10 

xyyx = -ry"^ - -rxy + -y+ - (4x^ - 3x - 1) . (101) 

Eqs. (78) and (79) yield in this case 



11 



-3z= (9.T + l)a;-3/5, (102) 

f(,) 6(x-l)(4. + l) 

^^'^^~ (6x-l)(9x+l)^- ^^""^^ 

Eq. (102) is a fifth order equation for x and seems to be transcendental. Obviously / (2) is not among the tabulated 
integrable functions. Eq. (101) does not coincide also with any of the equations in Sec. 1.3.4. from Ref. ||2^] which 
belong to the class (22). If we use Eq. (23) instead, g {() = l/f{z) and 

3Cx^/^ = 6x^ + 18x + l. (104) 

This equation is even more complicated than Eq. (102) and g (C) is not among the few integrable cases listed in 
Sec. 1.3. 2 from the same handbook. Unless some mistake has been made, Eq. (101) is integrable, but is not covered by 
Ref. P6[ . This situation is not unique. Recently, non-static charged perfect fluid distributions were discovered which 
too arc missing in the handbooks with solutions of differential equations |33] . 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have discussed explicit solutions of the Einstein equations for a static spherically symmetric perfect 
fluid with the 7-law equation of state. Three approaches may be found in the literature which, at first sight, have 
nothing in common. The approach of Collins transforms the TOV equation in polar gaussian coordinates into a 
two-dimensional autonomous system of differential equations. He studies it numerically H, |0], p9[ . This is a 
representative example of other similar dynamical systems [ |l6| , |lq ]. In his approach Klein finds a second-order 
autonomous master equation, whose solution determines all characteristics of the metric and the fluid. He lowers its 
order and studies it by series expansion and numerically. Klein also discovers a simple exact but singular solution. 
This work is done in curvature coordinates pffl , [0| . Haggag and Hajj-Boutros perform a similar study but in isotropic 
coordinates and only for a stiff fluid |14| . They also obtain an autonomous second-order equation and transform it 
into another second-order equation (see Eqs. (9) and (14) in their paper). Then they look for polynomial solutions 
and find either flat spacetime or the KT solution. 

We have generalized these three approaches and pursued them further, till their logical end - an equation, whose 
integrable cases are tabulated in the handbooks, like Ref. pM. Surprisingly, we invariably reach the Abel equation 
of the second kind (2). It behaves like a 'strange attractor' and underlines the common features in the different 
approaches. The question of integrability is answered by the form of / (z), the known integrable cases being tabulated. 
For this purpose we have generalized to an arbitrary coordinate system the dynamical system of Collins and to 
arbitrary n the HH approach . The functions z (x) and / (x) have the same general structure in all approaches but 
with different coefficients. The relations between the master variables x, xk and xhh have been elucidated, the 
TOV-CoUins approach serving as a basis. 

The integrable cases found are one and the same and include n = 00 (trivial dust solution), n = — 1 (de Sitter 
solution), n = —3 (Einstein static universe) and x = const, n < —5.83 or n > —0.17 (Klein- Tolman solution). 
They appear either as exceptional cases, when equations simpler than Eq. (2) are to be solved, or as Eq. (2) with 
/ (z) = az + (3, the simplest integrable case. Some candidate n lead to Eq. (63), but at least one of the constants 
there has the wrong value. 

It has been shown that Eq. (2) stands in the centre of the problem, independent from the approach or the coordinate 
system and the integrable cases may be derived in a unified manner. The problem is a very strange mixture of simple 
integrable cases and extremely difficult non-integrable ones. Perhaps this gives some explanation why the de Sitter 
and the Einstein solutions were found in 1917, two years after the appearance of general relativity, and why the next 
years have brought only a five-fold discovery of the KT solution. 

An important point is that the problem has been pushed to the mathematical realm of Abel differential equations 
and further progress depends on developments in this field. Using the Buchdahl transformation we have shown that 
the case n — —5 is integrable and have given its metric and fluid characteristics. However, its master equation in 
isotropic coordinates (101) does not seem to be present in the handbooks. This leaves this paper with an open end 
and one may hope that other integrable cases will be found in the future. 

Finally, it is interesting to note that all discrete integrable cases n = —1, —3, —5 are regular and fall in the interval 
where the singular KT solution does not exist. One is tempted to speculate that in this interval there is a regular 
one-parameter solution, encompassing the discrete cases. 
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